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We present a detailed calculation of the rare decay Z Hu'y'y via the quartic neutral gauge boson 
conpling ZZ 77 in the framework of the effective Lagrangian approach. The cnrrent experimental 
bound on this decay mode is then used to constrain the coefficients of this coupling. It is found that 
the bounds obtained in this way, of the order of 10 “^, are weaker than the ones obtained from the 
analysis of triple-boson production at the CERN collider LEP-2. 

PACS numbers: 12.60.Cn 


I. INTRODUCTION 

One of the most sensitive probes of physics beyond the standard model (SM) is provided by neutral gauge boson 
self couplings |Q. In the SM, trilinear and quartic neutral gauge boson couplings (TNGBCs and QNGBGs) ViVjVk 
and ViVjVkVi, with Vi = Z,7, vanish at tree level and their radiative corrections are known to be rather small, 
of the order of 10 “® - 10 “^° g, |j. Deviations from the SM predictions for TNGBCs and QNGBGs might point 
to new interactions such as those arising from stroimly-interacting electroweak models |^, a fourth family of chiral 
fermions with SM assignments of quantum numbers H , and the heavy fermions arising in the minimal supersymmetric 
standard model (MSSM) ||]. An interesting feature of QNGBGs involving at least one photon field stems from the 
fact that they are genuine in the sense that arise from effective operators that do not induce any trilinear gauge boson 
coupling. Therefore these genuine QNGBGs must be constrained from processes other than the ones used to constrain 
trilinear gauge boson couplings, such as boson-pair fusion or triple-boson production |^. This is to be contrasted 
with the case of the quartic lUIUyy coupling, which can arise from operators that also induce the trilinear WW^ 
coupling, which means that any constraint on the latter can be immediately translated into a bound on the former. 
Furthermore, while the quartic neutral ZZZZ coupling can be induced at tree-level, for instance by the exchange of a 
heavy scalar boson, QNGBGs involving at least one photon field can only arise at one-loop level or higher order in any 
renormalizable theory because of electromagnetic gauge invariance. In particular, QNGBGs have been constrained 
from the CERN e+e" collider LEP -2 data on triple gauge boson production [^. Considerable work has also been 
devoted to the analwis of QNGBGs through different processes at the next linear e+e” collider as well as 77, 67 and 
hadron colliders 

In this note we will present a detailed calculation of the ZZ77 coupling contribution to the rare decay Z —> vv^^. 
The current experimental limit on this decay mode 0J will then be used to constrain the coefficients of this QNGBC. 
To our knowledge, this calculation has never been presented in the literature. Our analysis will proceed in the same 
line as those presented in our previous works |^, where we obtained bounds on TNGBCs from the Z vv^ decay 
mode 1 ^ and on neutrino-photon interactions from Z —> [ p^ . We will find that the bounds obtained in this way 

rely on very few assumptions, though they are weaker than the ones obtained from the analysis of e“'"e“ —> .^77 and 
e+e" ^ W^W~^ data at LEP -2 i- The organization of the paper is as follows. In section II we present a short 
description of the effective Lagrangian for the ZZ77 coupling. Section III is devoted to present the calculation of the 
rare decay Z —> 1/1/77. Finally, our conclusions are presented in Section IV. 
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II. THE EFFECTIVE COUPLING ^^77 

When parametrizing physics beyond the Fermi scale in a model-independent manner by means of the effective 
Lagrangian technique there are two alternatives, i.e. the underlying new physics can be assumed to be either of 
decoupled or nondecoupled nature pd| ] . In the decoupling scenario it is assumed that the Higgs mechanism is realized 
in nature, thereby requiring the existence of at least one (relatively) light Higgs boson. In this case the low-energy 
theory is renormalizable d la Dyson, the decoupling theorem remains valid, and the effective Lagrangian is constructed 
out of operators respecting the SUl{‘ 2 ) x LV( 1 ) symmetry linearly. In this scenario the virtual heavy physics effects 
cannot affect dramatically the low-energy processes: the impact of new physics might become important only in 
those processes that are absent or very suppressed within the SM 0 - On the other hand, another possibility (the 
nondecoupling scenario) arises if the Higgs boson is very heavy or does not exist at all. There follows that the low- 
energy theory is nonrenormalizable due to the absence of the Higgs boson. This class of new physics can be assumed 
to be responsible for the symmetry breaking of the electroweak sector. In this scenario the SUl{‘ 2 ) x LV( 1 ) symmetry 
is nonlinearly realized[pT|. 

The lowest-dimension operators that induce the ZZ"f"f coupling have dimension six (eight) in the nonlinear (linear) 
realization. This means that any new physics effects arising from this coupling are likely to become more evident 
in the nonlinear scenario since in the linear one they are suppressed by higher powers of the new physics scale A. 
Therefore, in this work we will concentrate on the nonlinear scenario. Furthermore, only those operators that respect 
the custodial SUc{ 2 ) and the discrete C and P symmetries will be considered. It turns out that the operators that 
violate the custodial symmetry are tightly constrained by the p parameter. In the nonlinear scenario there are fourteen 
dimension-six operators that induce the ZZ^f^f coupling at tree level This was discussed to a large extent in 

Ref. and we will not dwell on this issue here. We rather focus on the Lorentz structure induced for the ZZ77 
coupling. In the unitarity ga.uge (U = 1 ), there are only two independent Lorentz structures for this coupling induced 
by dimension-six operators ||l2|] : 




16 A 2 , 


-aoF.,F^''Z^Z^ - 


-w 


16 A2 




Zr, 


( 1 ) 


where and = df_tZ^ — dyZ^. We have followed closely the notation introduced in |^. Below 

we will proceed to compute the decay Z —> izpyy using the above effective interaction. The experimental limit on this 
decay will then be used to constrain the coefficients oq/A^ and Cc/A^. 


III. THE RARE DECAY ^ ^ w'y'y 

We now turn to outline the calculation of the decay width for Z viyjj, which is somewhat similar to the one 
presented in Ref. for the rare decay Z vvAA in two Higgs doublet models, with A the CP-odd neutral 
scalar. In Ref. |l^ , the experimental bound on the rare Z —> 1/1/77 decay was used to constrain the neutrino-photon 
interactions and Dvjj. Here we will make an analysis along the same lines but focus on the purely bosonic ZZ77 
coupling, which contributes to the Z —> i/Djj decay through the Feynman diagram shown in Fig. 



FIG. 1: Contribution from the effective ZZyy coupling to the rare Z 1 / 1/77 decay. 
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The 4-vectors of the participating particles will be denoted as follows Z{p) —> v{pi)v{p 2 )Mki)^{k 2 ). The Feynman 


rule for the effective vertex Z{qfj_)Z{ql)^{ka)jik'^) is straightforwardly obtained from Eq 


2)liki 

■ i): 


g ^2 ^^^0 ik * k QaP k^kp) -\- dc k ' k T Pol0 (^/i k^ -\- ky k^ 

kp (^Pafi kjj T Pav kf^ k^ ippi^ kif ppy k^ 

where all the momenta are directed inward. The decay width can then be written as 


( 2 ) 


r(z ■ 


ni) = 


1 


(27r)®2^TO2 


iA4iv‘'(p-i:,.)n0. 

i=i i=i 


( 3 ) 


with Pi = pi, p 2 , ki and ^2 for i = 1, 2, 3 and 4, respectively. By the usual method and after a lengthy calculation 
we can obtain the squared amplitude. It reads 


\M\^ 


2 / _J_V ^t^i^PlP2 

3\Acw) ((p - fci - ^ 2)2 - m|)2 ’ 


( 4 ) 


where we have introduced the definition 

S;,!. = ^ (|4ao + acl^(fci • fc 2 )^ + 2 |acp (m| (fci • fe) + 2/ci • p ^2 • p) , (5) 

with Pfiv + 2p^py and = ki^k 2 i, + ki^k 2 fi- In Eq. a factor of 3 is included in the denominator as 

we are averaging over the Z boson polarizations; also, a factor of 2 has been introduced to account for two identical 
particles in the hnal state. The integration over pi and p 2 can be carried out straightforwardly with the aid of the 
following result Q 


1^''= /^^^^"^(Q-Pi-P 2 )p^P 2 " = J(QV" + 2 Q^Q"), ( 6 ) 

J Pi P 2 6 ' 

where Q = p — pi — P 2 - Once we are done with the integration over pi and p 2 , there still remains to integrate over ki 
and k 2 ■ To this end we will work in the center of mass frame of the Z boson. We hnd it useful to dehne the following 
variables ^ = 2p ■ = 2 k^jraz, ry = 2p ■ k2lm\ = 2 k^jmz, and w = (1 — cos0)/2. The decay width can thus 

be expressed as 


r(Z —> uv"f^) 




V A 2 , 


w 


mz a 


2 (2® 3 TT cw sw)^ 


h{^, rj, uj) dp duj, 


where 






(|4ao Tocp f{^,p,uj) + \ac\'^ p{^,p,uj)) , 


( 7 ) 


( 8 ) 


f{^,P,t^) = (12 + ^^ - p (12 - ry) - 2^ (6 - ry - 4a;p)) , 


(9) 


and 


p(^,p,a;) = 4 (1 + w) (1 + a; - 2a;f - 2pa; (1 — w^)). 


( 10 ) 
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It can be shown that the integration region is given by 11^ 


(11a) 

1-?7<C<1, ( 11 b) 

together with 0 < rj < 1. After numerical integration we are left with 

r(Z ^ l>i/ 77 ) = (lVoc|4ao + aeP + A^c|aeP) GeV, (12) 

with Nqc ~ 3.46 X 10“® and Nc ~ 10.31 x 10“®. In addition, if lepton universality is assumed, Eq. ( 0 ) is to be 
multiplied by 3 to account for all of the known neutrino species. From the LEP-2 data, the L3 collaboration set the 
following limit on .Z —> i /^77 


0 < w < 1 when 


?7 + ^- 1 


< w < 1 


when 


BR{Z Dvjj) < 3.1 X 10“®. (13) 

Assuming that either oq or Qc is dominant we obtain the following bounds 

oo » ttc, (14a) 

Qc » ao, (14b) 

which are weaker than those obtained at LEP-2 from ^77 and W^W~"f production ||]. However, it is important to 

note that the bounds based on the latter processes do not agree, as pointed out in Ref. Q, which means that much 
work along these lines is still required. In general Eqs. ( 0 ) and ( |l^ yield an allowed area in the oq /A^ versus ad A^ 
plane, as depicted in Fig. 


^ < 0.106 GeV^ if 
A2 

< 0.215 GeV^ if 


0.3 



-0.3 k . . , 

-0.15 -0.1 -0.05 0 0.05 0.1 0.15 
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FIG. 2: Allowed area (gray region) in the ao/A^ vs. a^/A^ plane from the experimental bound on the Z 1 / 1/77 decay. 

We can also take a different approach and instead of bounding the ZZ 77 coupling, we may use the most stringent 
bounds on it to predict its contribution to the Z —1 1 / 1/77 decay. From the most stringent experimental bound on 
the ao/A"^ and ad Id coefficients, of the order of 10 “^- 10 “^ §, it follows a limit on the contribution of the ZZ 77 
coupling to the Z i/i /77 decay: 


BR(Z ^ Z *77 ^ Pi/ 77 ) < 1 X 10"^^ (15) 

This indirect bound is above than the one found for the contribution of the neutrino-one-photon interaction Pi/ 7 , 
which is of the order of 10 “^^ p^ . 
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IV. CONCLUDING REMARKS 

In closing it is interesting to note that the ZZ77 coupling also contributes to the rare decay Z for which 

an experimental bound has already been set [p. This rare decay mode also receives contributions from the quartic 
^777 coupling. Therefore it would be possible, in principle, to use this rare decay to bound such QNGBCs. As far 
as the Z777 coupling is concerned, a tighter bound on it can be obtained from the three-body decay Z 777. The 
latter, together with the e+e" ^ 777 reaction, have been studied within the effective Lagrangian approach ||^ and 
we will not repeat the same analysis here. From the result presented in Ref. [T^ a bound on the effective Z777 vertex 
can be derived, which is of the same order of magnitude of the one obtained in this work for the ZZjj coupling. 

Finally we would like to emphasize that the importance of studying QNGBCs is rooted in the fact that they have 
a different origin than TNGBCs, i.e. they are induced by effective operators that do not induce any TNGBC. In this 
work we have found a constraint on the ZZjj coupling from the Z decay mode, which is weaker than the 

bounds derived from the analysis of ^77 and production at LEP- 2 . This is explained from the fact that 

the decay Z —> has a suppression factor due to the virtual Z boson propagator, which suffers less suppression 

when the initial Z boson is allowed to be off-shell. 
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